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We study a generalized Shastry-Sutherland model for the material SrCu2 (603)2- Along a line in the pa- 
rameter space, we can show rigorously that the model has a first order phase transition between Dimerized 
and Neel-ordered ground states. Furthermore, when a magnetic field is applied in the Dimerized phase, 
magnetization plateaus develop at commensurate values of the magnetization. We also discuss various as- 
pects of the phase diagram and properties of this model away from this exactly soluble line, which include 
gap-closing continuous transitions between Dimerized and magnetically ordered phases. 



In recent years, many novel magnetic materials have 
been synthesized which exhibit spin-gap behavior. In 
these materials the ground state is a spin-singlet and 
there is a gap to all spin excitations. Such phenom- 
ena have long been studied in quasi-one dimensional sys- 
tems, but much recent interest has arisen from the discov- 
ery of quasi-two dimensional spin-gap materials CaV409 
§, NazTizSbsO || and SrCu2(B03)2 @- The latter 
material is particularly interesting in that, by virtue of 
the crystal geometry, it is an experimental realization 
of the Shastry-Sutherland model B, for which an ex- 
act dimerized singlet eigenstate can be written down, 
which for a range of parameters is the ground state of 
the model. Among the interesting experimental findings 
for SrCu2 (603)2 are that the system appears very close 
to a transition to a Neel phase and it also shows magne- 
tization plateaus as a function of magnetic field [^|| . 




Here we consider a generalized Shastry-Sutherland 
model, with Hamiltonian 



H — Ji 2_^ ^i ' ^j + "^2 2_^ Si ■ Sk + J3 2_^ Si ■ Si. 
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FIG. 1. The geometry of spins and couplings in the mate- 
rial SrCu2(B03)2. 



where the bonds corresponding to interactions Ji , J2 and 
J3 are shown in Fig. 1. We assume Ji is antiferromag- 
netic and can henceforth set Ji to unity. For J3 = 
the model reduces to the original Shastry-Sutherland 
model. For J2 = J3, the model has infinitely many 
conserved quantities. The total spin on each Ji bond 
commutes with TL and thus each eigenstate of the Hamil- 
tonian can be characterized by the number and position 
of the triplets present. These triplets then form (in gen- 
eral a diluted) spin-one Heisenberg model, with nearest- 
neighbor interactions on the square-lattice. It is easy to 
show that this model has three phases at T = with 
first order transitions between them. For large negative 
J2 the ground state is a fully polarized ferromagnet, for 
large positive J2 the ground state is equivalent to the 
Neel ground state of a spin-one Heisenberg model on the 
square-lattice, which is rigorously known to be ordered 
H, and whose numerical properties are very well estab- 
lished 0]. In between the ground state is the Dimerized 
Shastry-Sutherland singlet phase. At the Dimcr to Neel 
transition the sublattice magnetization jumps from to 
about 80 percent of the classical value. 

Away from the J2 — J3 line, we use symmetry argu- 
ments, Dimer series expansions M together with consid- 
erations of the classical limit to discuss the ground-state 
phase diagram and properties of this model. The model 
with J2 and J3 interchanged can be mapped into the 
original one by interchanging the spins on all Ji bonds. 
Thus the phase diagram is symmetric with respect to 
the J2 = J3 line. We will concentrate our discussion on 



the region J2 > J3. First, let us compare energies of 
various classically ordered phases with the energy of the 
Dimer phase to get a first view of various phases and 
their boundaries. 
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FIG. 2. Phase diagram for the model showing ferromag- 
netic (F), Neel (N), columnar (Cf and Cjv), helical (H) and 
dimer (D) phases. See text for discussion of phase boundaries. 

The different phases are shown in Fig. 2, with the 
Dimer phase denoted by D. The magnetically ordered 
phases are best described with respect to a rotated lat- 
tice where individual spins have four nearest neighbors 
(connected by J2 bonds) |J]. This lattice is topologi- 
cally equivalent to a square-lattice. On this lattice, in 
addition to a ferromagnetic phase (F) and an antiferro- 
magnetic Neel phase (N), there are columnar (Cat and 
Cf) and helical (H) phases. The two columnar phases 
are colinear and equivalent to each other. In the Cjq 
phase the spins order antiferromagnetically along one of 
the axes and have period four (ttii) in the perpendicular 
direction. In the Cp phase, the spins are ferromagnetic 
in one direction and have period four in the other direc- 
tion. Classically the four helical phases can be mapped 
onto each other. In the helical phase for J2 > — J3 > 0, 
the helix runs along one of the axes, with successive spins 
rotating by an amount +9 as one moves along that direc- 
tion H . Along the perpendicular direction the change in 
spin directions alternates between -\-9 and —9. The angle 
9 is non-unique and is one of the solutions to the equation 



cos (tt — 9) 



2(J2 - J3 



Jl + ^Ji - 24J3(J2 - J3) 



(2) 



Defining x — J2/J1, y — J^/Ji, the columnar, helical, 
Dimer triple point in the classical phase diagram is lo- 
cated at xtr = (45 - 4V6)/36 « 0.9778, ytr = -(9 -I- 
4V6)/36 « -0.5222. The Neel-helical phase boundary is 



given by X -|- 5?/ = 1, whereas the asymptotic {x — > 00) 
phase boundary between the helical and columnar phases 
is given by the equation y = -^(— 8a; -f 1 -f OiXjx)). 

Since there are no quantum fluctuations in the ferro- 
magnetic and Dimer phases and on the boundary be- 
tween ferromagnetic and helical phases, these remain 
true phase boundaries in the model and are shown by 
solid lines in Fig. 2. Note that the Dimer- ferromagnetic 
boundary is first order, whereas the Dimer-helical bound- 
ary is second order. The other classical phase bound- 
aries are shown by dashed lines. They leave an oval-like 
Dimer phase in the middle. The first-order Dimer-Neel 
phase boundary can be determined quite accurately along 
Ji — Jz to be at x = 0.42957(2), and has been given in 
previous numerical studies along J3 = (and equiva- 
lently J2 = 0) g. These points are shown by solid dots 
and we connect them smoothly to indicate a first order 
Dimcr-Neel phase boundary. 

Actually, it is not evident whether the Dimer-Neel 
phase boundary along the line J3 = is first or second 
order, or even whether there is an intermediate phase be- 
tween the two. Albrecht and Mila ||] have argued that 
there is an intermediate helical phase between the Neel 
and Dimer phases. Their Schwinger Boson mean- field 
treatment leads to the estimate that the Neel phase ex- 
tends only down to x « 0.91, whereas the helical phase 
exists between 0.61 < a; < 0.91. On the other hand, 
using series expansions, Weihong et al. H have argued 
that the Neel phase extends down to x = 0.691 at which 
point there is a first order transition to the Dimer phase. 
The finite-size calculations fipl also do not suggest any 
helical phases. Though, Albrecht and Mila have argued 
that this is because the helical phases are not properly 
accomodated in finite geometries. 

The quantitative validity of Schwinger Boson calcula- 
tions is difficult to judge. One generally expects quantum 
fluctuations to stabilize colinear phases. And this could 
considerably reduce the extent of the helical phases in the 
phase diagram. In several spin models, where numerical 
calculations have been done, the sublattice magnetiza- 
tion of the Neel phase goes continuously to zero and it 
is separated from incommensurate phases by a singlet 
phase ||ll|]. On general grounds, Ferrer |1J] has argued 
that the Neel phase must be separated from helicoidal 
phases by an intermediate spin-liquid phase. Thus, it is 
reasonable to assume that along J3 = there is a direct 
transition between Dimer and Neel phases. 

Along J3 = 0, Weihong et al. estimate that the Dimer 
to Neel transition happens at J2/J1 = 0.691(6). Using 
d-log Fade approximants to analyze the gap series, we 
estimate that it vanishes at J2/J1 ~ 0.697(2). Thus, 
within the uncertainties of the series analysis, this tran- 
sition could be continuous. Around this value of the cou- 
plings, the sublattice magnetization series from the Neel 
side is also consistent with zero [pi. The primary rea- 
son for believing that the transition is first order is that 



the energies for the Neel and Dimer phases appear to 
cross at a non-zero angle. However, if the transition is 
continuous, then very close to the transition, the Neel 
energy curves should change slope p. Thus, given all 
the numerical evidence, a plausible conclusion is that the 
transition is very weakly first order, though it could also 
be second order. 

Given the above, it is natural to expect some continu- 
ous transitions when J3 < 0. To explore this possibility, 
we have developed series expansions for the triplet ex- 
citations in the dimer phase to 15th order for arbitrary 
Jz/J2 using the flow equation method ||I^. For J3 — 0, 
the expansion coefficients agree with those calculated by 
Weihong et al. The unit cell of the lattice contains two 
dimers, giving rise to two triplet modes. These modes 
are almost degenerate throughout the Brillouin zone and 
become exactly degenerate as g — > due to symmetry. 
We find that the gap minimum is always at q equal to 
zero even as one moves from the Neel towards the fer- 
romagnetic phase. We use d-log Fade approximants to 
calculate the locus of points, where the triplet gap closes. 
This contour is also depicted by a solid line in Fig. 2. It 
marks a boundary at which the Dimer phase becomes 
locally unstable, and hence the dimer phase cannot ex- 
ist beyond that line. Without studying all eigenstates of 
our system, it is not possible to say if some other level 
crossing transition leads to a different ground state before 
we get to this line. It is plausible to think that at least 
parts of this line represent continuous phase transitions 
between the Dimer and magnetically ordered phases. As 
seen in Fig. 2, a possible consistent scenario is that very 
near the J3 = line, we have a multicritical point where 
a second order transition line meets a first order phase 
boundary. 

These continuous phase transitions between the Dimer 
phase and various magnetic phases are rather unusual. 
They are not in the conventional 2 -I- 1-dimensional 0(3) 
universality class as expected for the non-linear sigma 
models |l3|. This is evident from the fact that in the 
Dimer phase, the ground state remains unchanged and 
hence the correlation length remains of order unity. In 
contrast, for a generic dimerized spin system, the correla- 
tion length gradually grows and diverges as the gap goes 
to zero u%- The continuous phase transition, here, is 
somewhat analogous to the density driven generic phase 
transitions in the Bosonic Mott insulators [|15| . 

However, there are some important differences. Un- 
like the case of Bosonic Mott insulators, the spectrum 
appears to be linear at the transition. Along the J3 = 
line, we estimate that the gap vanishes at a; = 0.697(2), 
with an apparent exponent v of 0.45(2). Different d- 
log Fade approximants show remarkable consistency with 
each other. Fig. 3 shows the spectrum, in the reduced 
Brillouin zone, at the transition. Different ways of ana- 
lyzing the series all point to a finite spin-wave velocity 
and a linear spectrum. These results suggest that this 



transition belongs to a new universality class W] 
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FIG. 3. The triplet spectrum along the J3 = line at 
X = 0.697, where the gap vanishes. 

With our present calculations we cannot study the 
transitions from the ordered side and thus cannot estab- 
lish the full nature of these phase transitions nor can we 
say anything about the stability of columnar and helical 
phases in the overall phase diagram. Quantum ffuctua- 
tions can lead to additional singlet (spin-gap) phases be- 
tween the Neel and the helical phases and possibly elim- 
inate helical phases altogether from the Neel side of the 
phase diagram. These Neel to singlet phase transitions 
should be similar to those found in the Ji — J2 square- 
lattice Heisenberg model @ . 
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FIG. 4. Magnetization as a function of magnetic field. Note 
that the plateaus are valid for J4 < J2 < 0.43Ji. 

When a magnetic field is applied to the Dimer ground 
state along the special line J2=J3, the resulting magneti- 
zation is shown in Fig. 4. The triplet excitations, aligned 
by the field, have no dispersion, but a nearest neighbor 
repulsion. Thus they form a simple Wigner crystal (or a 
Bosonic Mott insulator) at one half the saturation mag- 
netization. If we add an additional weak antiferromag- 



netic coupling between neighboring horizontal Ji bonds 
(and similarly neighboring vertical Ji bonds) of the form 
JiiSi + S2) ■ {S3 + S4), the triplets remain dispersionless 
but they now have an additional second neighbor repul- 
sion. This leads to additional Wigner crystal phases at 
one and three quarter fillings and hence plateaus in the 
magnetization at one and three quarters of the satura- 
tion value. In Fig. 4, we also show the ordering pattern 
of the Wigner crystal on different plateaus by showing 
four Ji bonds. A line denotes a {S^ — 1) triplet on the 
bond whereas a circle denotes a singlet. 

As we move away from the J2 = J3 line, the triplet 
develops dispersion. It is useful to think of the prob- 
lem in terms of the Bose Hubbard model [n6[, with the 
(S"^ = 1) triplet representing hard core Bosons. These 
Bosons have a strong nearest-neighbor repulsion and a 
weak diagonal (second-neighbor) and further-range hop- 
ping. At half-filling, the strong repulsion will clearly lead 
to a Wigner crystal and hence the magnetization plateaus 
will remain. However, the transition between the magne- 
tization plateaus may now be second order. Indications 
of such plateaus also exist in the finite size calculations of 
Miyahara and Ueda for J3 = Q . However, due to the 
finite size, all plateaus appear discontinuously in their 
study. The question of whether there will be additional 
magnetization plateaus at other rational fillings perhaps 
including valence bond states, as in one dimension p9[ , 
deserves further attention. 

We now comment on the materials: One would naively 
expect the SrCu2(B03)2 system to be close to J3 = 0, a 
limit that has been considered by other authors ^M . The 
ratio of J2 to Ji has been placed in the literature ^ close 
to the Dimer to Neel transition. Even though this transi- 
tion may be first order, it would be very weakly so, due to 
the vicinity of the multicritical point. In this sense, this 
material may allow us to study a special quantum criti- 
cal point, not the generic Neel to singlet quantum phase 
transition. However, this transition may be unstable to 
the generic 0(3) transition, when the special eigenstate 
of Shastry and Sutherland is not a true eigenstate due 
to some small perturbations. It would be interesting to 
study this crossover theoretically. An interesting prob- 
lem could be the instability of these special transitions 
due to spin-lattice couplings. 

Magnetization plateaus have been observed in the ma- 
terial at one eighth and one quarter of the saturation 
magnetization. The exactly soluble model suggests that 
these phases may be regarded as simple Wigner crystals 
of local triplets. The question of whether the models 
away from J2 = J3 will have magnetization plateaus at 
other rational fractions, or whether other couplings such 
as J4 are needed for this, deserves further attention. 
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